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Bulk motion measurements in clusters of galaxies with ATHENA-like 
missions 
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The hierarchical formation of clusters of galaxies by accretion of material releases gravitational energy which dissipates 
into the intracluster gas. The process heats the material and generates gas turbulence and bulk motions and thus kinetic 
pressure. Mapping the velocity fields of the moving subunits would enable a new diagnostics tool for cluster formation 
studies and unbiased X-ray mass estimates. The required spatially resolved high resolution spectroscopy is not currently 
available. I demonstrate here the feasibility of detecting and mapping the velocities of the bulk motions using the Doppler 
shift of the Fe XXV Ka line with the proposed ATHENA satellite. 
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1 Introduction 

During its life cycle, a cluster of galaxies experiences col- 
lisions and mergers with other clusters and subunits. Even- 
tually the relaxation processes take place and the intraclus- 
ter material approaches hydrostatic and virial equilibrium. 
Consequently, there will be bulk motions in the cluster ma- 
terial at different distance scales and velocity levels, de- 
pending on the magnitude of the event and at which phase 
we observe the cluster 

Mapping the velocity fields of the bulk motions would 
open a new tool for studying the dynamics of the clus- 
ters. Comparing the velocity maps with those predicted by 
cosmological simulations would be useful for the study of 
the formation of the large scale structure and provide con- 
straints on cosmological parameters. Combining with the 
measured galaxy density field, the bulk velocity maps could 
be used to test whether the fundamental picture of gravita- 
tional collapse is coiTect (Branchini et al. 2001; da Costa et 
al. 1998; Dore et al., 2003). 

The bulk motions produce kinetic pressure which can 
reach a level of 10% level of that of the thermal gas even in 
most relaxed clusters (Lau et al. 2009). If not accounted for, 
the kinetic pressure may bias the hydrostatic X-ray mass es- 
timates low by ~10% at rsoo- Thus, it would be important 
to measure this component in clusters in order to obtain un- 
biased mass estimates for the cosmological applications. 

In case of recent, strong collisions which happen close 
to the plane of the sky, the merger shocks can be observed by 
the X-ray moiphology, as in the rare cases of A520 (Marke- 
vitch et al. 2005), A754 (Macario et al. 2011) and A2146 
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(Russell et al. 2012). However, in most lines of sight the 
merger features are hidden by the projection. This allows 
the possibility of measuring the radial bulk motions via the 
Doppler shift of the emission lines. 

Currently the constrains on radial bulk motions are 
rather poor due to the lack of spatially resolved high spec- 
tral resolution instruments. The proposed ATHENA mis- 
sion would have carried such an instrument, an X-ray Mi- 
crocaloiimeter Spectrometer XMS. In this paper I will ex- 
amine the expected quality of bulk motion measurements 
with a mission approximating the capabilities of the pro- 
posed ATHENA instruments, i.e. an ATHENA-like mission. 

2 Bulk motions in clusters 
2.1 Range of velocities 

It is widely accepted that clusters of galaxies form by merg- 
ing of smaller structures of matter accreted along large scale 
structure filaments. The velocities of the accretion flows are 
assumed to reach a level of 1000 km s~^ (e.g. Frenk et al. 
1999). A collision of two clusters or protoclusters of simi- 
lar mass produces a strong merger The related bulk motion 
velocities can reach several 1000 km s^^ as in the case of 
the Bullet cluster (e.g. Markevitch et al. 2002). Simulations 
of Nagai et al. (2002;2003) indicate that ongoing minor 
mergers may produce bulk velocities up to ^1000 km s^^. 
The released gravitational energy is dissipated into the in- 
tracluster material which eventually approaches the hydro- 
static equilibrium. However, simulations indicate that even 
in the most relaxed clusters there are residual bulk motions 
present throughout the cluster volume at 100 km s~^ level 
(e.g. Lau et al., 2009; Nagai et al. 2003). 
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2.2 Doppler shift 

The radial components of the bulk motions of 100 - 1000 
km s~^ correspond to 2-20 eV shift of the Fe XXV and 
XXVI Ka emission line centroid energies (~6 keV). It is 
challenging to use the currently most powerful X-ray instru- 
ments at 6 keV (XMM-Newton/EPIC, Chandra/ACIS and 
SUZAKU/XIS CCDs) for these measurements due to lim- 
itations in the energy resolution (^100 eV). The Gaussian 
centroid can still be determined with precision better than 
100 eV, assuming that the instrument gain is accurately cal- 
ibrated. This requirement can be relaxed if one considers 
relatime motions between the main cluster and the moving 
part. In the following I discuss only the statistical precision 
of the line centroid determination, and not the calibration 
issues nor the uncertainties of the cosmic redshift measure- 
ments. 

2.3 Observational constraints 

Suzaku/XIS instruments have been used to place upper lim- 
its for the bulk motion velocities at ^1000 km s^^ level in 
several clusters: A2319 (Sugawara et al. 2009), Centaurus 
(Ota et al., 2007), AWM7 (Sato et al. 2008). A first sig- 
nificant detection of a bulk motion has been achieved with 
Suzaku for the merging subclump in A2256 (Tamura et al. 

201 1) . The radial velocity difference between the main clus- 
ter and the subclump is 1500±300±300km s^^ (where the 
two sets of uncertainties refer to statistical and systematical 
ones, respectively). 

3 ATHENA-like missions 

ATHENA (Advanced Telescope for High ENergy Astro- 
physics) was one of three L-class (large) missions being 
considered by the European Space Agency in the Cosmic 
Vision 2015-2025 plan. In May 2012 the Jupiter mission 
Jupiter Icy Moons Explorer (JUICE) was chosen for launch. 
However, ESA has committed to continue supporting tech- 
nology developments for a future large X-ray facility. At the 
time of writing this paper (Oct 2012) there was an under- 
standing that ESA would shortly appoint a small team from 
the community to provide input based on the the ATHENA 
study team activities. 

In this paper I used the responses and background esti- 
mates for the X-ray Microcalorimeter spectrometer (XMS) 
and a wide field imager (WFI) as reported in the ATHENA 
assessment report used by the ESA Space Science Advi- 
sory Committee (the ATHENA Yellow Book, Barcons et al., 

2012) . 

3.1 X-ray Microcalorimeter Spectrometer XMS 

The requirement for the energy resolution of XMS is 3 eV 
at 6 keV. With its 100-1000 times higher effective area at 
0.5 keV, compared to those of the current high resolution 



instruments onboard XMM-Newton and Chandra, and spa- 
tial resolution of 10 arcsec, XMS would enable spatially 
resolved high resolution spectroscopy. This would yield a 
breakthrough in mapping the spectral properties for ex- 
tended sources like clusters. Also, the bandpass extends to 
12 keV which allows the measurement of the Fe XXV K a 
emission line, prominent in clusters of galaxies. The down- 
side of XMS is the relatively small FOV (2.3x2.3 arcmin^) 
which renders the velocity mapping of a whole cluster chal- 
lenging. 

3.2 Wide Field Imager WFI 

Even though the energy resolution of WFI (150 eV at 6 keV) 
is not sufficient to resolve the components of Fe Ka com- 
plex, it can still determine very precisely the centroid of the 
Gaussian distribution of the total line emission (see Fig. [U 
due to high photon statistics. This is due to the very large 
effective area (^0.5 m^) of the X-ray telescopes onboard 
ATHENA. Thus also WFI is a powerful tool for mapping 
the cluster velocities. 



T ' ' ' ' 1 ' ' ' ' 1 ' ' ^ 

A2256 lOOks WFI simulation 



Fe XXV Ka centroid = 6.311 keV 
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Energy [keV] 



Fig. 1 Data from a lOOks WFI simulation of the A2256 
main cluster (red crosses) and the subclump (blue crosses) 
together with the input models (solid curves). The verti- 
cal dashed lines highlight the centroids of the Fe XXV Ka 
lines. 



4 Improvement with ATHENA-like missions 

I examine here the capability of the proposed ATHENA 
instruments for measuring the expected bulk motions 
in clusters. I used the MEKAL model (Kaastra 1992) 
in the XSPEC package to model the cluster emission 
(bremsstrahlung continuum + coUisional excitation lines) 
adopting the metal abundances of Grevesse & Sauval 
(1998). 

I used the current estimates of the instrument responses 
(Barcons et al., 2012) to simulate the cluster spectra, includ- 
ing the instrument background and 90% resolved cosmic X- 
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Table 1 Velocity precision for a cluster with z = 0. 1 and 
kT = 5 keV with lOOks observation 
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ray background. Then I fitted the simulated data in the 5.5- 
7.5 keV band (excluding the channels where the Fe XXVI 
Ka emission is significant) with a model consisting of a 
power-law component for the continuum and a Gaussian 
line for the Fe XXV Ka emission. The 1 a statistical un- 
certainty of the Gaussian centroid then yielded the estimate 
for the statistical precision of the velocity measurement, as 
summarised in Table [T] 

4.1 A2256 and the subclump 

In order to obtain a real-life example of the performance 
of XMS and WFI for measuring the velocities of a nearby 
minor merger, I used the Suzaku results for the A2256 and 
the subclump (Tamura et al. 201 1) to simulate spectra using 
an exposure time of 100 ks. 

A fit to the simulated XMS data (see Fig. ^ yielded 
a statistical uncertainty of the redshift (cr^ '--^10^^) corre- 
sponding to a velocity precision at a level of 1 km s^^, 
i.e. a 1500(7 detection for the clump motion. Measurement 
is very precise because many line features are resolved and 
each centroid gives weight to x^. Dividing the XMS emis- 
sion of the subclump into boxes with width of 30 arcsec 
(i.e. 5x5 pixel map for the full XMS FOV) yielded a sta- 
tistical precision level of ^^10 km s^^. This level of detail 
would provide a breakthrough in modelling the dynamics of 
the merging subclumps. Since the flux of A2256 subclump 
is comparable to that of the central region of similar size in 
A2256 (Tamura et al., 201 1), the above calculations yield an 
approximate estimate of the expected velocity mapping pre- 
cision level in the bright nearby cluster centres with XMS. 

A spectral fit to the data from a 100 ks WFI simulation 
of A2256 main cluster and subclump (see Fig. [U yielded a 
statistical precision for Fe XXV Ka shift corresponding to 
velocity precision of 60 km s^^, i.e. a detection at ^25 cr 
level. 

4.2 Mapping the whole cluster 

In order to derive general conclusions about the perfor- 
mance of ATHENA-like missions for the velocity mapping 
in clusters, one should simulate spectra with a grid of repre- 
sentative values for cluster temperatures, metal abundances 
and redshifts and realistic exposure times. Also, one should 
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consider a range of luminosities for the moving regions and 
their velocities and directions of motion. 

In this work I limited the complexity of the above ap- 
proach by using an exposure time of lOOks for a cluster 
with kT = 5 keV and z=0.1, assuming that all of the emis- 
sion in the line of sight of a studied region originates from 
the moving subunit. I further assumed that the luminosity 
of the moving subunit at a given distance from the cluster 
centre can be estimated with a a /3 - model for the surface 
brightness with (3 = 2/3 and rcoro = 0. 1 rsoo. This may be an 
underestimate of the actual signal since the likely enhance- 
ment of the emission due to the subunit is not accounted for. 
I adopted a bolometric luminosity Lf,o;(i"50o) = 7 x 10'*'* erg 
s^^ within rsoo from L-T relation of Pratt et al. (2009). Us- 
ing rsoo - T relation of Vikhlinin et al. (2006) I adopted rsoo 
= 10 arcmin. 



ioorA2256 subclump lOOks XMS simulation- 




6.30 6.32 6.34 6.36 

Energy [keV] 



Fig. 2 Data from a lOOks XMS simulation of the A2256 
subclump (blue crosses) together with the input model, 
shifted by 100 km s^^ (red line). 



4.2.1 XMS 

A single XMS pointing covers only a fraction of the full 
cluster volume: at z = 0. 1 the rsoo radius for a cluster with 
kT = 5 keV is ^1 Mpc (Vikhlinin et al., 2006) and covers 
an area of ^300 arcmin^, i.e. ~^50 higher than the FOV of 
XMS. Thus it is not feasible to cover the full r = rsoo re- 
gion in nearby clusters with XMS. Limiting the mapping 
into the central r=0.25 rsoo region would require ^ 10 point- 
ings which might be feasible for a few clusters. Single XMS 
pointings at larger radii might still be useful, and I thus cal- 
culate estimates for these observations in the following. 

I simulated XMS spectra with the above generic clus- 
ter (see Section|42]i within the fuH FOV i.e. (0.2 rsoo)^- 
Fitting the simulated spectra yielded a statistical precision 
of -200 km s-i at i^0.25 rsoo- At r=0.5 rsoo ^90% of the 
total signal in the 5.5-7.5 keV band is due to background 
emission which degrades the statistical precision of veloc- 
ity to —400 km s^^ and at rsoo the data are so noisy due to 
the background that reliable velocity measurements cannot 
be obtained. 
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If the above cluster was located at z=l .0, a single central 
XMS pointing would cover the cluster out to ^0.5 r5oo. The 
received flux would be similar as that of the z=0. 1 cluster at 
0.5 rsoo, when using the full XMS FOV. Thus, with a central 
lOOks XMS pointing it is feasible to obtain a single bulk 
motion measurement at z = 1.0 assuming the full cluster 
within 0.5 r5oo is moving with a radial velocity component 
higher than 400 km s~^ . 

4.2.2 WFI 

Fitting the simulated data using the generic cluster de- 
scribed in Section 14.21 showed that in the centre the bulk 
velocity can be mapped with angular resolution of '-^O.l 
rsoo with statistical precision of ^ 100 km s^^. At r=0.25 
rsoo the lower cluster flux requires a larger extraction re- 
gion to achieve similar statistical quality as in the centre. 
Using an extraction box with size of (^0.2 rsoo)^ yields a 
velocity precision of ^ 200 km s^^. At a distance of 0.5 rsoo 
from the centre, the background dominates (as in the case of 
XMS) and the velocity precision degrades to ~800 km s^^ 
level when using an extraction box size of ^0.3 rsoo)- Re- 
peating the exercise with a very hot cluster (kT =10 keV) 
improves the continuum signal but the Fe XXV emission 
decreases and thus the velocity precision does not improve 
significantly. 

5 Conclusions and discussion 

The simulations showed that for a nearby (z < 0.1) cluster 
with kT = 5 keV, using a 100 ks exposure using the instru- 
ments proposed for ATHENA, one could 

- measure bulk motions in cluster centres at ~ 0.1 rsoo 
spatial scale with statistical velocity precision level of 
-10 (100) km s-i using XMS (WFI) 

- map cluster velocities in ~ 10 regions with WFI with 
a single central pointing up to r = 0.25 rsoo with —200 
km s^^ statistical precision with angular resolution de- 
creasing from — 0.1 rsoo to — 0.2 rsoo with increasing 
radius 

- obtain —200 km s^^ precision with XMS up to 0.25 
rsoo, but to cover the fuU cluster within this radius one 
needs — 10 pointings 

- obtain a velocity precision of — 400 km s^^ with a sin- 
gle XSM off-axis pointing at 0.5 rsoo using the emission 
from the fuU FOV 

- obtain a velocity precision of — 800 km s^^ with WFI 
at 0.5 rsoo with a single central pointing 

- not obtain meaningful velocity measurements at rsoo 
due to dominating background 

The above estimates for the level of statistical precision 
of velocity measurements with XMS and WFI indicate that 
an ATHENA-like mission would enable for the first time the 
mapping of radial components of the bulk motions due to 
recent minor and major mergers in nearby (z<0.1) clusters 



of galaxies within the central 0.25 rsoo- Also the residual 
motions due to past mergers can be mapped in these cen- 
tral regions. Single velocity measurements at the distance 
of 0.5 rsoo for a cluster at z=0. 1 at the level of —400 km s~^ 
can be achieved with lOOks off-axis pointings using XMS. 
At higher distances the mapping is more limited due to the 
background but the average motion within the central 0.5 
rsoo up to z=l .0 can be measured with a central lOOks point- 
ing with XMS, if the radial velocity component is bigger 
than -400 km s"^ 

These measurements would open a new tool to study 
cluster dynamics. This would be complementary to the cur- 
rent analyses of intracluster shocks and sloshing and thus 
give a better handle on the 3-dimensional dynamics of clus- 
ters. This would yield a breakthrough in modelling the clus- 
ter physics and its implications to cosmology. 

The kinetic Sunyaev Zeldovic (kSZ) effect as measured 
by WMAP (KashHnsky et al., 2008; Osborne et al., 2012) 
has been used to constrain bulk motions of clusters. How- 
ever these results have been controversial due to relatively 
weak kSZ signal compared to the thermal one. Thus, one 
must combine a large sample of clusters and assume a sim- 
ilar motion for all the clusters in the sample. Mak et al. 
(2011) calculate that a sample of 400 clusters would al- 
low Planck to detect a coherent 500 km s^^ bulk motion. 
Clearly ATHENA would make improvement here, being 
able to measure the velocities in individual clusters with bet- 
ter precision and spatial resolution. 
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